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We present a theoretical model for the experimentally found but counter-intuitive exceptionally 
long lifetime of surface nanobubbles. We can explain why, under normal experimental conditions, 
surface nanobubbles are stable for many hours or even up to days rather than the expected mi- 
croseconds. The limited gas diffusion through the water in the far field, the cooperative effect of 
nanobubble clusters, and the pinned contact line of the nanobubbles lead to the slow dissolution 
rate. 



Introduction - Since their first prediction and discov- 
ery almost 20 years ago intense research on surface 
nanobubbles has raised many questions about this in- 
triguing and important phenomenon which has great po- 
tential for various applications Surface nanobub- 
bles have now been widely reported on various surfaces 
in contact with water employing various detection mech- 
anisms like Atomic Force Microscopy (AFM) and most 
recently also through direct optical visualization 0, 0). 
With all these different methods they are found to be- 
have differently than regular macroscopic bubbles. Sur- 
face nanobubbles behave peculiar in several ways: their 
contact angle is always much lower than expected from 
Young's law [1, Q , they are stable against violent decom- 
pression (lo| . and in particular they are stable for much 
longer than expected: For such small bubbles one would 
expect a lifetime of order /xs, due to the high Laplace 
pressure inside the bubbles which drives the gas into the 
liquid. On this last question many explanations were pro- 
posed, ranging from contamination that shields or limits 
the diffusive outflux of gas to a dynamic equilibrium 
situation where lost gas is replenished IH, 13 1. However, 
both theories are refuted by experimantal evidence: the 
addition of surfactants does not influence the behaviour 
of nanobubbles [3] , and the circulatory gas flow required 
for the dynamic equilibrium theory is not measured in 
all experiments so it cannot be the stabilization mecha- 
nism 0, HH . In addition, a large problem with the dy- 
namic equilibrium theory is that it requires some form of 
driving to satisfy the second law of thermodynamics, and 
its origin is unclear. In molecular dynamics, some local 
inflow near the contact line was indeed observed, but its 
strength was too weak to explain the stability of surface 
nanobubbles 
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A different approach is therefore required, in this Let- 
ter we provide an alternative explanation for the long 
lifetimes of surface nanobubbles. The theory relies only 
on classical, well-known, and proven continuum concepts 
such as diffusion and Henry's law. Furthermore, the the- 
ory only uses confirmed nanobubble properties, namely 
the pinned contact line [ly, LLZf of nanobubbles and the 
fact that nanobubbles exist in high coverage fractions at 




FIG. 1. (Color online) Sketch of a liquid layer of thickness 
I in contact with a solid (left). The top of the liquid is ex- 
posed to atmospheric conditions. At the solid-liquid interface 
nanobubbles are present with typical internal contact angle 
6, height h, and radius of curvature R. They are not drawn 
to scale. The arrows indicate the gas flow direction. On the 
right a further enlargement of one nanobubble is shown. 



the liquid-solid interface 0, [l9|. No fitting or uncon- 
trolled assumptions are required to obtain lifetimes that 
are consistent with experimental findings. 

The Letter is organized as follows. First, the theory is 
explained and the relevant equations are derived. Next, 
we solve the equations numerically and analytically. We 
vary several parameters to demonstrate the robustness 
of the long lifetimes of surface nanobubbles in varying 
experimental conditions. In addition, we apply the the- 
ory to the case of electrolytically generated nanobubbles 
and find that also here it is consistent with experimental 
results. We conclude with predictions from the theory 
which can straightforwardly be tested in experiments. 

Theory - We consider an infinitely large plate in con- 
tact with a liquid layer with thickness I (Fig. [IJ. The 
liquid layer is in contact with the atmosphere at z = £ 
and the solid-liquid interface is located at z = 0. The 
solid-liquid interface is covered with nanobubbles with a 
number density (per area) of p, hence the (average) spac- 
ing between neighbouring nanobubbles is (c?b„b) ~ -^=. In 
experimental studies nanobubbles are always recovered in 
high coverage densities [18l42l| . Assuming quasi-steady 
diffusion, any variation of the dissolved gas concentration 
4> in the horizontal (x, y) direction due to non-uniform gas 
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outflux decays as exp(—irz/db u b)- Hence wc can assume 
that for £ > 5 (db u b) the the diffusion of gas through the 
liquid layer is governed by the one-dimensional diffusion 
equation: 



at 



D 



dz 



z e [0,£] 
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with D the diffusion constant of the gas in the ambient 
liquid and <p the (number) concentration of gas in the 
liquid. The characteristic timescale for diffusion through 
the water layer is £ 2 /D ~ f0 5 s, which is similar to the 
lifetimes obtained in experiments. The vast difference (ixs 
vs. days) compared to previous estimates originates in 
using £ as the relevant length scale instead of the bubble 
radius R, i.e. we use the far-field length-scale. Using R as 
the length-scale is justified for a (free) bubble in an infi- 
nite medium where. However, in the case of nanobubbles 
it is important to realize that gas has not left the system 
until it is released into the atmosphere and hence £ is the 
relevant length scale as this is the distance the gas has 
to travel through the liquid from the bubble towards the 
atmosphere. This is also apparant in a thought experi- 
ment where we do not allow gas to leave the liquid into 
the atmosphere, for example by putting the liquid drop 
in a closed container. Since the liquid is supersaturated 
with some gas that has left the bubble, and which cannot 
escape, an equilibrium is reached (described by Henry's 
law discussed later in this work) and the bubbles do not 
dissolve |22| . When opening the container this excess dis- 
solved gas can be released into the ambient air, allowing 
the bubble to lose gas into the liquid. This 'traffic jam 
effect' plays a crucial role also in the long nanobubblc 
lifetimes. 

We will now fully describe the boundary conditions 
for above differential equation (TTJ). They are given by 
Henry's law, which relates the gas concentration in the 
liquid to the gas pressure outside the liquid near the in- 
terfaces: 

Hz = o, t ) = Pb " b(J ? (jV(f))) and <f>(z = £, t) Patm 
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Here, ku is Henry's constant, Pbub the pressure inside the 
nanobubbles and p a tm the atmospheric pressure. R(t) is 
the radius of curvature of the nanobubbles, which de- 
pends on time because the bubbles get flatter as they 
drain. In the case of a pinned contact line, the radius 
of curvature is related to the (internal) contact angle 
0(t) by R(t) = si ^g(t) , where Rb is the base radius (cf. 
Fig. [TJ which is constant due to the pinned contact line. 
The (relative) pressure inside the bubbles Pbubit) is the 
Laplace pressure 
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Pbub{t) = Patm + -5- = Patm + 77" sill 9(t) , (3) 
it rib 

where 7 is the liquid- vapour surface tension. For 8 < 90° , 
which is the case for surface nanobubbles, the internal 



pressure thus decreases as 9 decreases. This effect pro- 
vides a negative feedback in the dissolution process, pro- 
longing the lifetime of the nanobubbles. In this work we 
do not consider the effects of electrostatic effects on the 
internal pressure of surface nanobubbles, as electrostatic 
effects act to reduce the internal pressure and are there- 
fore not a driving force but rather a stabilizing force. It is 
therefore possible that the derived lifetimes in this work 
(hours, days) are an underestimation of real lifetimes. 

The (single) nanobubble gas content N(t) decreases 
due to the diffusive flux of gas at z = 0, which is the 
location at which gas is transfered from the nanobubbles 
(z = _ ) to the liquid (z = + ). The diffusive flux is 
given by Fick's law J = —D(d<fi/dz), thus 
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The factor p arises to convert the molecular flux per unit 
area of the substrate to the molecular flux per single sur- 
face nanobubble. Eq. ((4]) immediately shows how a low 
nanobubble coverage {p ~ 0) corresponds to small global 
flux, J = -pdN/dt. 

To evaluate the boundary condition at z = (Eq. |J2])), 
we need to relate the number of atoms inside a nanobub- 
ble (Eq. (jU)) to the geometrical shape of the bubble. 
To calculate the geometrical properties of a nanobub- 
ble containing N atoms of gas, we use the ideal gas law 
p(0)V(8) = NksT using the expression for the volume of 
a spherical cap V{9) = § (i? b /sin6») 3 -(2-3cos(9+cos 3 6>). 
Here, kg is Boltzmann's constant, T the temperature (as- 
sumed to be constant at 300K), and 9 the gas-side (in- 
ternal) contact angle of the nanobubbles. This implicit 
equation can be solved numerically for 9(N). We can 
then calculate the radius of curvature of the nanobubbles 
R, which gives us the internal bubble pressure [Eq. fl^J)]. 
Using this pressure, the boundary condition at z = 
(Eq. ([3])) can be evaluated, which closes our model. In 
the next section, we will solve these model equations nu- 
merically. 

Numerical evaluation - Due to the non-trivial bound- 
ary condition at z = (Eq. ©) we first solve the diffusion 
eq. (fTJ) numerically. The simulations were done for differ- 
ent initial conditions. Since the real initial conditions are 
unknown, we choose the two extremes between which we 
expect the real initial conditions. The first type of ini- 
tial conditions consists of a linear concentration profile, 
which allows the system to begin transporting gas from 
the bubbles immediately (t = 0), 
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Here, Rq indicates the initial radius of curvature of the 
nanobubbles. We choose Ro such that it is equivalent to 
an initial contact angle 9q = 40° for given base radius 
Rb- 
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FIG. 2. (Color online) Results of the calculations using the 
initial conditions from Eq. ([5| and parameters described in 
the text, (a): Snapshots of the concentration profile at 1 
hour intervals. Since the bubbles are pinned, the radius of 
the curvature increases as the bubble drains, lowering the 
concentration at the bubble side (z — 0). (b): Evolution of the 
nanobubble gas content (black) and the radius of curvature 
of the liquid-gas interface R (red) as a function of time. 



The second type of initial conditions assumes that 
the nanobubble formation procedure (ethanol-water ex- 



change or replacing cold water with warm water 21[) 
supersaturates and mixes the water such that the con- 
centration is uniform and equal to the concentration near 
the nanobubbles: 
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The real initial concentration profile will most likely be 
something between ((SJ) and ([5]): the ethanol-water ex- 
change uniformly supersaturates the water but it takes 
the nanobubbles some time to form so some gas already 
drains into the atmosphere. As we will see, both initial 
concentration profiles produce long-living nanobubbles. 

Results - How long docs a nanobubble survive ac- 
cording to this description? Using typical parameters 
that apply to experiments on surface nanobubbles Q 
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D = 10" 9 m 2 /s, k H = 2.6 • 10~ 19 Pa-m, and T = 300 K) 
and the initial condition ([5]) we obtain the results shown 
in Fig. [2] Fig. [2^ shows hourly snapshots of the con- 
centration profile (j)(z). From these curves it is apparent 
that the transport of gas away from the bubble is limited 
by the diffusion rate of the gas through water far away 
from the bubble. This leads to an almost flat (zero-slope) 
concentration profile near the bubble through which the 
diffusive gas flux is very small. Fig. (red) shows the 
radius of curvature of the spherical cap R as a function 
of time. This radius of curvature increases (due to the 
pinned contact line) which lowers the internal gas pres- 
sure [Eq. <j3j) ] , enhancing the lifetime of the nanobubbles. 
Finally, Fig. (black) shows the amount of particles 
contained inside an individual nanobubble through time. 
Here we see that a typical nanobubble which starts out 
with over 100,000 atoms and ends up with just over 500 



atoms after 36 hours. 

As a criterion for bubble dissolution we choose a crit- 
ical bubble height of h* = 1 nm, after which it may 
no longer be appropriate to use continuum physics, and 
also effects such as the disjoining pressure start to domi- 
nate the pressure inside the bubble rather than solely the 
Laplace pressure. From molecular dynamics it is known 
that continuum models (e.g. Navier-Stokes, diffusion 
equation, Henry's law) are valid down to the nanometer 
scale 15J, |22( . Using the same criterion we find that above 
bubble is stable for 36 hours, which is 10 orders of mag- 
nitude longer than previously thought and in agreement 
with experimental findings and the expected timescale 

i 2 /d. 

Analytic solution — The concentration profiles in 
Fig. [2] suggest that the boundary condition at z = 
is of Neumann-type, d<f>/ dz\ z= a = 0. Taking the time 
derivative of the boundary condition at z = [Eq. ((2)] 
and substituting Eq. (j4]) into the result we obtain 
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tion is 4>{z,t) = (/) bubfi e^p(—K 2 Dt/(4£ 2 )) 
^transient (*)■ The exact form of A 



Filling in representative values for the quantities we 
find that the concentration gradient at z — is over 
5 orders of magnitude smaller than the typical global 
concentration gradient (« <f>/£), with (f> = <fi — (f) a tm- 
This means that the boundary condition for the gra- 
dient at z = can indeed be considered to be ap- 
proximately zero. For zero gradient the analytic solu- 

COs(tTz/(2£)) + 

transient depends on the 
initial conditions, but declines quickly. A remarkable fea- 
ture of this result is that <f>{z,t) only depends on D and 
£, and is completely independent of p and 7. Of course, 
p must be high enough to be able to consider the system 
as 1-dimcnsional. Similarly, the fraction p/ [Dd<j)/ dN) 
must be low enough such that the local gradient at z = 
[Eq. 0] is small compared to the global gradient. 

Varying initial conditions How does the initial con- 
centration profile affect the lifetime of the nanobubbles? 
During the ethanol-water exchange procedure that is 
most commonly used to generate nanobubbles experi- 
mentally, ethanol is flushed away with clean water. It is 
therefore likely that the initial gas concentration profile 
is uniform in z, due to mixing. We redid the same calcu- 
lations as before, using an uniform initial concentration 
profile, the results are plotted in Fig. |3l We observe very 
similar concentration profiles as before, except for small 
times where the influence of the initial conditions is still 
felt. As can be observed in Fig. [3] (b), it takes 2-3 hours 
before the bubbles 'feel' the influence of the ambient air 
and start to dissolve. This means that these first hours, 
the bubbles barely shrink as the (global) concentration 
gradient near the bubbles is close to zero. 

Robustness of the results - How robust arc the 
nanobubble lifetimes against variations in the experimen- 
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FIG. 3. (Color online) Numerical results for the initial condi- 
tions JfjJ. (a): Hourly snapshots of the concentration profile. 
Due to the initial conditions, some gas that is initially dis- 
solved in the liquid first has to drain to the atmosphere until 
a concentration gradient exists near the bubble. As compared 
to the linear profile case (Fig. [2J the bubbles gain an addi- 
tional few hours of lifetime because of this, (b): Evolution 
of the nanobubble gas content and the radius of curvature of 
the surface nanobubbles. 



tal system? By changing 7 and I in the numerical calcu- 
lation, we verify the result from the analytical solution to 
the diffusion equation for which it holds that the diffusion 
profile evolution only depends on D and I. 

First, we look at the surface tension 7. Most of- 
ten an ethanol-watcr exchange procedure is applied to 
form nanobubbles. This method introduces contamina- 
tion into the system which lowers the surface tension. It 
is therefore important to understand the influence of 7 
on the nanobubble lifetime. We find that surface ten- 
sion does not play any role in the dissolution time of 
nanobubbles. This result remains counter-intuitive as 
surface tension is the driving force for nanobubble dis- 
solution. Indeed, a higher surface tension increases the 
Laplace pressure [Eq. (J3J], thus increasing the driving 
that leads to dissolution. However, it also increases the 
gas content inside the bubbles. This denser reservoir re- 
quires a larger flux to drain in the same time. Both effects 
scale linearly with 7, hence they cancel. 

In previous studies, the liquid height £ has never been 
considered as a parameter that affects nanobubble life- 
times. Based on the analytical solution presented before 
we expect that r ~ £ 2 /D, which is indeed exactly re- 
covered from the numerical calculations. This indicates 
that the liquid cell size (and geometry) is very important 
to the lifetime of nanobubbles: the amount of liquid the 
gas has to travel through determines the timescale of the 
nanobubble lifetime. £ is easy to vary, we suggest to per- 
form corresponding experiments to test our prediction. 

Electrolysis - How do nanobubbles behave according 
to this description when gas is generated at the solid- 
liquid interface (e.g. by electrolysis 3, 23l)? A rough 
estimate based on the values cited in ref. [l6| gives a con- 
stant influx of order ~ 10 5 molecules per bubble per sec- 



ond. For an equilibrium to exist (the bubbles neither 
grow or shrink) the diffusive flux away from the bubble 
must then be equal to this influx due to electrolysis, hence 
—<j>(z = 0)/£ = fO 5 s _1 . This corresponds to a nanobub- 
ble contact angle of 9 = 49°, and nanobubble content 
N = 2 • fO 5 . Interestingly, this means that the bubble 
contents are refreshed every two seconds. This only high- 
lights the fact that nanobubbles are not static: without 
driving (such as electrolysis) they dissolve whereas with 
driving the gas atoms inside the bubble are replaced ev- 
ery couple of seconds. 

Conclusion - We have modelled the gas flow from 
nanobubbles through the liquid to the atmosphere to 
study the lifetime of surface nanobubbles. We find that 
nanobubbles are not stable, but dissolve by diffusion. 
However, due to their pinned contact lines and that the 
gas has to diffuse towards the atmosphere, they dissolve 
on a much longer timescale than free bubbles in an in- 
finite liquid. This last point (diffusion through a rela- 
tively thick layer of liquid) explains why previous molec- 
ular dynamics results could not recover long nanobubble 
lifetimes: the system size in molecular dynamics is lim- 
ited to several tens of nanometers, resulting in lifetimes 
of order 100 ns, consistent with our findings in ref. [15. 
Using the correct length-scales (£ ~ 1 — 10 mm), we find 
that surface nanobubbles can easily survive in excess of a 
day, an increase of 10 orders of magnitudes as compared 
to previous estimates of fis. The results from the exper- 
iment where electrolysis is used corroborate the results 
in this Letter: a non-zero current is measured meaning 
that gas is continuously formed at the substrate. We 
find that the gas flux induced by this current leads to 
nanobubbles with 9 = 46°, consistent with experimental 
results of ref. 



1G 



J.H. Snoeijcr, H. Gcldcrblom, and X.H. Zhang are grate- 
fully acknowledged for discussions. This work is part of 
the research programme of the Foundation for Funda- 
mental Research on Matter (FOM), which is part of the 
Netherlands Organisation for Scientific Research (NWO). 



[1] J. Parker, P. Claesson, and P. Attard. Bubbles, cavities, 
and the long-ranged attraction between hydrophobic sur- 
faces, J. Phys. Chem. 98, 8468 (1994). 

[2] M. A. Hampton and A. V. Nguyen, Nanobubbles and the 
nanobubble bridging capillary force, Adv. Colloid Inter- 
face Sci. 154, 30 (2010). 

[3] J. R. T. Seddon and D. Lohse, Nanobubbles and mi- 
cropancakes: gaseous domains on immersed substrates, 
J. Phys.: Condens. Matter 23, (2011). 

[4] V. S. J. Craig, Very small bubbles at surfaces-the 
nanobubble puzzle, Soft Matter 7, 40 (2011). 

[5] J. R. T. Seddon, D. Lohse, W. A. Ducker, and V. S. J. 
Craig, A Deliberation on Nanobubbles at Surfaces and in 
Bulk, ChemPhysChem 13, 2179 (2012). 



5 



[6] S. Karpitschka, E. Dietrich, J. R. T. Seddon, H. J. W. 

Zandvliet, D. Lohse, and H. Riegler, Nonintrusive Optical 

Visualization of Surface Nanobubbles, Phys. Rev. Lett. 

109, 066102 (2012). 
[7] C. U. Chan and C.-D. Ohl, TIRF Microscopy for the 

Study of Nanobubble Dynamics, Phys. Rev. Lett., in print 

(2012). 

[8] B. M. Borkent, S. de Beer, F. Mugele, and D. Lohse, 
On the Shape of Surface Nanobubbles, Langmuir 26, 260 
(2010). 

[9] X. H. Zhang, N. Maeda, and V. S. J. Craig, Physical 
Properties of Nanobubbles on Hydrophobic Surfaces in 
Water and Aqueous Solutions, Langmuir 22, 5025 (2006). 

[10] B. M. Borkent, S. M. Dammer, H. Schonherr, G. J. 
Vancso, and D. Lohse, Super stability of Surface Nanobub- 
bles, Phys. Rev. Lett. 98, 204502 (2007). 

[11] W. A. Ducker, Contact Angle and Stability of Interfacial 
Nanobubbles, Langmuir 25, 8907 (2009). 

[12] M. P. Brenner and D. Lohse, Dynamic Equilibrium Mech- 
anism for Surface Nanobubble Stabilization, Phys. Rev. 
Lett. 101, 214505 (2008). 

[13] J. R. T. Seddon, H. J. W. Zandvliet, and D. Lohse, Knud- 
sen Gas Provides Nanobubble Stability, Phys. Rev. Lett. 
107, (2011). 

[14] X. Zhang, M. H. Uddin, H. Yang, G. Toikka, W. Ducker, 
and N. Maeda, Effects of Surfactants on the Formation 
and the Stability of Interfacial Nanobubbles, Langmuir 
28, 10471 (2012). 

[15] J. H. Weijs, J. H. Snoeijer, and D. Lohse, Formation 
of Surface Nanobubbles and the Universality of Their 
Contact Angles: A Molecular Dynamics Approach, Phys. 



Rev. Lett. 108, (2012). 

[16] S. Yang, P. Tsai, E. S. Kooij, A. Prosperetti, H. J. W. 
Zandvliet, and D. Lohse, Electrolytically Generated 
Nanobubbles on Highly Orientated Pyrolytic Graphite 
Surfaces, Langmuir 25, 1466 (2009). 

[17] X. Zhang, D. Chan, D. Wang, and N. Maeda, Stability 
of Interfacial Nanobubbles, Langmuir, submitted (2012). 

[18] S. Yang, S. M. Dammer, N. Bremond, H. J. W. Zandvliet, 
E. S. Kooij, and D. Lohse, Characterization of Nanobub- 
bles on Hydrophobic Surfaces in Water, Langmuir 23, 
7072 (2007). 

[19] S. Yang, E. S. Kooij, B. Poelsema, D. Lohse, and H. J. W. 
Zandvliet, Correlation between geometry and nanobubble 
distribution on HOPG surface, Europhys. Lett. 81, 64006 
(2008). 

[20] J. W. G. Tyrrell and P. Attard, Atomic Force Micro- 
scope Images of Nanobubbles on a Hydrophobic Surface 
and Corresponding ForceSeparation Data, Langmuir 18, 
160 (2002). 

[21] X. H. Zhang, X. Zhang, J. Sun, Z. Zhang, G. Li, H. Fang, 
X. Xiao, X. Zeng, and J. Hu, Detection of Novel Gaseous 
States at the Highly Oriented Pyrolytic GraphiteWater 
Interface, Langmuir 23, 1778 (2007). 

[22] J. H. Weijs, J. R. T. Seddon, and D. Lohse, Dif- 
fusive Shielding Stabilizes Bulk Nanobubble Clusters, 
ChemPhysChem 13, 2197 (2012). 

[23] L. Zhang, Y. Zhang, X. Zhang, Z. Li, G. Shen, M. Ye, 
C. Fan, H. Fang, and J. Hu, Electrochemically Controlled 
Formation and Growth of Hydrogen Nanobubbles, Lang- 
muir 22, 8109 (2006). 



